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Abstract 

We report on the perturbative 0(as) correction to the light-cone QCD sum rule for 
the i? —> TT transition form factor /"*■. The correction to the product /s/'*' in leading 
twist approximation is found to be about 30%, that is similar in magnitude to the 
corresponding 0{as) correction in the two-point sum rule for fs- The resulting 
cancellation of large QCD corrections in /■*“ eliminates one important uncertainty 
in the sum-rule prediction for this form factor. 
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1. The semileptonic decay B ^ ttIui is one of the most important reactions for the 
determination of the CKM parameter Vub- However, in order to extract Vub from data 
one needs an accurate theoretical calculation of the hadronic matrix element 

{7i{q)\uj^b\B{p + q)) = 2/+(/)g^ + (/+(/) + (1) 

where p + q, q and p denote the B and vr four-momenta and the momentum transfer, 
respectively, and are two independent form factors. 

A very reliable approach to calculate in the framework of QCD is provided by the 
operator product expansion (OPE) on the light-cone BBi in combination with QCD 
sum rule techniques. The sum rule for the form factor /^(p^) has been obtained in [^, 
taking into account all twist 2, 3 and 4 operators, while is derived in 0. The most 

important missing elements of these calculations are the perturbative QCD corrections 
to the correlation function leading to (i- Here we report on a calculation of the 0(a<j) 
correction to /■*“ which eliminates one of the main uncertainties in the existing sum rule 
results. 

The calculation has several aspects which are worth pointing out. Firstly, the sum rule 
is actually derived for the product /b/^, fs being the B meson decay constant dehned 
by 

(H|6i75d|0) = m%fB/mb . (2) 

The form factor /+ itself is then obtained by dividing out fs taking the value determined 
from the corresponding two-point QCD sum rule. In previous estimates, the 0{as) correc¬ 
tion to Jb was thereby ignored for consistency because of the lack of the 0(a<j) correction 
to /b/^- Our calculation now allows to take into account the correction to /b which 
is known to be sizeable. Secondly, knowing the 0{as) corrections, also the heavy quark 
mass entering the sum rule can be properly dehned. Thirdly, perturbative corrections to 
exclusive amplitudes involving light-cone wave functions have so far been studied only 
for massless quarks. For example, in ^ ^ the amplitude of the pion transition to two 
virtual photons was calculated to 0{as)- The calculation for a hnite quark mass is new 
and will have numerous applications. 

The main result of our work is the following. The 0{as) correction to the light-cone 
sum rule for the product /b/’*' calculated in the leading twist approximation is about 30% 
and positive. Since the 0{as) correction to /b is similar in size and of the same sign, the 
large QCD corrections cancel in f~^ making the prediction of the form factor very reliable, 
at least from the point of view of perturbative QCD. 

In this letter, we outline our calculation, present the hnal analytical results, and give 
hrst numerical estimates. Technical details, a thorough numerical analysis, and further 
applications will be presented elsewhere. 

2. The light-cone sum rule for the form factor f~^ was derived in [Q. We repeat here 
the necessary points. In QCD, the correlation function of two heavy-light currents, 

Ft,{p,q) = i j dxe^P'^{Tr{q)\T{u{x)'yf,b{x),mbb{0)i'y5d{0)}\0) 

= {P + (P + q?)Pf^ , ( 3 ) 
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can be calculated in the region {p + < 0 and < ml — O(lGeV^) using OPE near 

the light-cone, i.e. at ~ 0. In (^, we have multiplied the pseudoscalar current by 
the 6-quark mass in order to assure renormalization-group invariance of the correlation 
function. After contracting the 6-quark helds in (^, is expressed in terms of bilocal 
matrix elements of increasing twist. In the present calculation we focus on the leading 
twist 2 contribution which enters through the following matrix element: 


7r(g) 


M(x) 7 ^ 75 Pexp iwsj da X ■ A{ax) J (i(0) 


j + ■■■ , (4) 

0 


where the ellipses stand for terms of higher twist. The path-ordered gluon operator ensures 
gauge invariance. In the light-cone gauge, a; • A = 0, adopted here as usual this operator 
is unity. The distribution function (pn{u) is known as the twist 2 light-cone wave function 
of the pion 0, I, i. 

Comparison of (|^ and (|) shows that in order to calculate /■*■ one only has to deal 
with the invariant amplitude P in (|^). With (^, P can be written as a convolution of 
a hard amplitude T(p^, {p + qYiU) calculable within perturbation theory, with the pion 
wave function ipT,{u) containing the long-distance effects: 

F{p^,{p +qf) = -U j dup^{u)T{p^,{p + qf,u). (5) 

JO 

In zeroth order in the hard amplitude represented graphically in Fig. la reads 


^o(p^(p + g)^M) 


__ 

p^(l — u) + {p + qYu — ml 


( 6 ) 


At hxed p^ < ml, F is an analytic function in the complex [p -|- g)^-plane, with a cut 
along the real axis starting from {p + q)'^ = ml- 
One can therefore write a dispersion relation 


F{pMp + qf) 


r°o s)ds 

ml s-{p + g)2 


(7) 


Equating the QCD result obtained with = ^ImP and the hadronic representation 

of P following from (|^) with the spectral density 

s) = 6(s - m|)2m|/s/+(/) + s)0(s - sq) (8) 


yields the desired relation between f~^ and the invariant function P. In (|[), the hrst term 
stems from the B ground state, whereas the second term represents the contributions from 
the higher resonances and the continuum in the P-meson channel above the threshold sq. 
Invoking quark-hadron duality the latter is replaced by the spectral density . The 
sum rule hnally follows from the above after Borel transformation in {p + qY'- 




1 

2m\ 



QCD 


(p^s)e 



ds , 


(9) 


3 








where 


pQC'D(p 2 ^^) ^ f s, m) . (10) 

71 J 
0 

With the zeroth order approximation (^, one easily obtains 

2 

ImTo(p^, s, u) = — 7 r 5 (l — -^(1 — u) - %u). (11) 

mi mi 

Substitution of (|T 0 |) and (pd]) in (^ and integration over s reproduce the leading twist 2 
contribution to the light-cone sum rule given in Q. In this approximation the evolution 
of is taken into account in the leading order (LO). In order to go to the next-to-leading 
order (NLO), one has to calculate the 0{as) correction to ImT and use the NLO-evolution 
of (Ptt- This problem is solved below. 


3. The hrst step is to calculate the 0{as) correction to the hard amplitude T which 
we write as 

ol C 

T{ri,r2,u) = To{ri,r2,u) + ^^^^Ti(ri,ra,u) , (12) 

dvr 

introducing convenient dimensionless variables ri = and r 2 = {p + qY/m^. The 

zeroth order amplitude Tq is given in (P). In Figs, lb - Ig we show the Feynman dia¬ 
grams determining the hrst order amplitude Ti. The calculation is performed in general 
covariant gauge in order to have a possibility to check the gauge invariance of the re¬ 
sult. Both the ultraviolet (UV) and infrared divergences are regularized by dimensional 
regularization and renormalized in the MS scheme with totally anticommuting 75 . This 
choice is motivated by the fact that the same scheme is used in the calculation of the 
NLO evolution kernel of the wave function (Pt,{u) |T0[| . 

From the diagrams depicted in Fig. 1 we hnd 


Ti(ri,r2,u) = 


with 


+ 


3(1+ P) 
(1-p)^ 
2 

(ri - r2Y 


mt 


A-ln^ + ll- 

, /i^ ) 1 - p L 

1 — / X ^ / X 


u 


G{p)-G (n) 


2G(p)-G(ri)-G(r2) 
1 — ri r 


1 — M L 


G{p)-G (r 2 ) 


p-h (1 - p) ln(l - p) 2 (1 -r 2 )ln(l-r 2 ) ^ 3-p 


p^ 


1 -p 


r2 


Y-pY 


( 1 -p)ln(l-p) ( 1 -r 2 )ln(l-r 2 )' 


(1 - M)(ri - ra) 


P 


r2 


(13) 


A = 


4-d 


- 7 e-F ln( 47 r), p = n + u{r 2 - n), 


mz 


G{p) = Li 2 (p) + ln (l-p)-ln(l-p) ( A-ln—+ II , 


( 14 ) 
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IA 2 {x) = — / y ln(l —t) being the Spence function. The UV renormalization scale and the 
0 

factorization scale of the collinear (COL) divergences are taken to be equal and denoted 
by /i. In order to trace the origin of the various divergent terms we have performed 
additional explicit calculations. In particular, we have used mass regularization by giving 
the light quarks a small but hnite mass, and momentum regularization keeping the light 
quarks off mass shell. In this way, we have unambiguously separated the COL-divergent 
terms from the UV-divergent terms. The latter add up to 

rr(>'i.r2.«) = (15) 

The correlation function (^) involves the unrenormalized quark currents J 5 = and 
Jfj_ = U'jfib as well as the bare 6 -quark mass m^. As usual, we dehne the corresponding 
renormalized quantities by 


J5 ^ -^ 5>^5 5 Ju ^ ^vJui ^ ^mdlb 


In the MS-scheme, the renormalization constants are given by 


^5 


1 + 3A 


asCp 

Att 


Zv = 1 , Zn 


1 - 3A 


asCp 

An 


(16) 


We see that the overall renormalization factor of is Z^Z^Zy = 1. The UV-renormalized 
hard amplitude T then follows from the unrenormalized result (|^) just by reexpressing 
the unrenormalized mass rrih through the renormalized mass mi, . As a result, an ad¬ 
ditional 0{as) contribution to T emerges which exactly cancels the term given in 


In addition to the UV-divergent terms, the function Ti contains the COL-divergent 
terms: 


T^^^{ryr2,u) = -AToiu) 


3 - 21n 

-2 In 


1 - r 2 \ (1 - ri)(l - r 2 ) - u{l - ri)(r 2 - ri) 


1 — ri 


u{r2 - riy 


(17) 


1 - p (1 - ri)(l - r 2 ) - u{l - u){r 2 - rif 
1 - r2 


u{l — u){r 2 — riY 


It is straightforward to check that can be written in the form 

1 U 


T^^^{ri,r2,u) = -A- [ dwV{w,u)To{ri,r2,w) 

2 Jo 


(18) 


where V{w, u) is the kernel of the Brodsky-Lepage evolution equation @ of the light-cone 
wave function ipT,{u) introduced in ( 0 ): 


di~p,,{u, ^)/dha.^ = / dijjV{u,ijj)^p,,{ijj, 

Jo 


with 


V{w, u) = 




TT 


1 — W 
. 1 — M 


1 + 


w — u 


w 


Q{w — u) -\— 1 - 1 - 


u 


u — w 


0(m — w) 


(19) 


,( 20 ) 
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The operation + is defined by 


1 

V{w,u)+ = V{w,u) — 5{w — u) Jv{v,u)dv. 


( 21 ) 


The appearance of COL-divergent terms in the hard amplitude T in the form ([18|) reflects 
the factorization of the correlation function into a wave function and a hard amplitude 
Bii- For the factorization scheme we have again adopted the MF-scheme, i.e. we have 
subtracted the terms in the UV-renormalized hard scattering amplitude proportional to 
A. These are the terms absorbed in the definition of the scale-dependent wave function. 
The remaining /i-dependences of the hard scattering amplitude and of the wave function 
compensate each other. 

Up to now we have worked in the MF-scheme. However, the MS quark mass depends 
explicitly on the renormalization scale fi and implicitly on the renormalization prescrip¬ 
tion. A renormalization-scheme-independent definition of the quark mass within QCD 
perturbation theory is given by the pole mass which we denote ml. Since we intend to 
use the set of parameters {ml, fs, Sq) determined self-consistently from an independent 
analysis of the two-point sum rule for fs [|^ it is convenient to replace by ml also in 
the sum rule for /b/"*" using the well-known one-loop relation: 


rhb = ml [l + 


Air 


-4 -|- 3 In 


m 


*2 ' 




( 22 ) 


To 0{as), this replacement adds the term 


2p 


a-py 


4 — 3 In 


mf' 




(23) 


to the renormalized amplitude Ti. The final result for the invariant function (j^) then 
reads 


i 

F{ri.r2) ^ ~f^ J du<p„{u. fi)T(ri.r2.u, . 


(24) 


where the renormalized hard amplitude is given by 


T(ri,r2,M,p) = 



[2G (p) — G {ri) — G (r2)] 


/ X2 V (p) ~ (’^i)] + “j- - [G {p) - G {r2)] 

(^’i — ^ 2 ) \ u 1 — u 

p-h (1 - p) ln(l - p) 2 (1-r 2 )ln(l-r 2 ) 2 


P" 


+ 


p- 1 


r2 


p- 1 


'(1-p)ln(l-p) (1-r 2 )ln(l-rs)' 


(1 - u){ri - r 2 ) 


P 


r2 


(25) 
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Here G{p) = G(p)|a =05 and ip.„{u, p) is the pion wave function evolved to the scale p in 
NLO. 

To proceed further according to (J^) and (p!0D we calculate the imaginary part of the 
hard scattering amplitude (P^D for r 2 > 1 and ri < 1 : 

1 


-Imr(ri,r 2 ,M,/x) = <5(1 - p) 

71 


as{p)CF 


dvr 


+ 2 Li 2 (l - r 2 ) -2 (In 


r2 


5{l- p) 

1.2 


*2 

TT^ - 6 + 31n^^^^^^^-2Li2(ri) 


p^ 


1 — ri 


+ 2 In r 2 + 


r2 


r2 


+e{p-l) 


8 


ln(p - 1 ) 


P 


+ 21 nr 2 H- 2 — 2 ln(r 2 — 1 ) + hi 


(2 ln(r 2 — 1 ) — ln(l — ri)) 
mt^\ 1 


r2 


p^ 


P 


-2 


r2 


{n - r2){p - ri) 


In p — 2 ln(p — 1 ) + 1 — In 


m 


*2 ^ 


p 


+2 


1 — ri 


(ri - r 2 )(r 2 - p) V r 2 

+ 0(1 - p) 


In 2 ^ - 2 In ^^ 


r2 


y lap 

p- 1 


0 1 A 

1 ^ 

1 

1 ] 

+ 2 1 


H —2 

— 

r2- p\p 

r2) 


P\ 


/ JJl 

2 In r 2 H- 2 ln(r 2 — 1 ) — In 


*2 ' 


r2 


p^ 


P 


ri 


(ri - r 2 )(r 2 - p) 

Here, the operation + is dehned by 

1 


In r 2 + 1 — 2 ln(r 2 — 1 ) — In 


m 


*2 ' 


r2 


P 


r 2 - p r 2 


(26) 


dpfip) 


P 


= J dp{f{p) - f{l)) 


P 


(27) 


This prescription takes care of the spurious infrared divergencies which one encounters by 
taking the imaginary part of (| 2 ^) . 

Substituting (^BD and (PDD to (j^) one obtains the desired sum rule in 0{as) for the 
form factor f~^ in the leading twist 2 approximation: 






277m\ 


rsQ 


ds J du (p^(M,p)ImT(-£^, 


2 


'ds 


(28) 


The subleading twist 3 and 4 contributions are presently known only in zeroth order in 
as They will be taken into account in the numerical analysis. 


4. The second step is to determine the decay constant fs and the pion wave function 
p) in NLO. For that purpose we have analyzed the two-point sum rule for fs 
obtained from the renormalization-group-invariant correlation function 
ml{0 I T{J^(a;)+ 5 ( 0 )} | 0) in 0 (q;s) [|^]. For the running coupling constant we use the 
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two-loop expression with Nf = 4 and = 234 MeV corresponding to as{Mz) = 


0.112. For /i^ we take the valne /i^ = — ml^ corresponding to the average virtnality 

of the correlation function which in turn is given by the Borel mass parameter M^. With 
this choice the following correlated results are extracted from the two-point sum rule: 


fe = 180 d= 30 MeV 


m. 


= 4.7^0.! GeV, Sq = 35 ± 2 GeV^ 


(29) 


In the following, we adopt the central values in the above intervals. Note that without 
0{as) correction one obtains Jb = 140 ± 30 MeV. The remaining parameters entering 
(p 8 |) are directly measured: rriB = 5.279 GeV and = 132 MeV. 

The wave function ip-„ can be expanded in terms of Gegenbauer polynomials 
G))/^( 2 n — 1). Arguments based on conformal spin expansion allows one to neglect 


higher terms in this expansion. We adopt the ansatz suggested in ||I41] : 

/To) = 4'o(m) a2(/io)^2(M) + a4(/io)^4(M), 


(30) 


where T„(m) = 6 n(l —m)G^/^(2m —1). The asymptotic wave function ipn{u) = 6u{l — u) 
is unambigously hxed 0. The terms n > 0 describe nonasymptotic corrections. The 
coefficients 02 (/xo) = 2/3 and 04 (po) = 0.43 at the scale /xo = 500 MeV have been extracted 
from a two-point QGD sum rule for the moments of (pn{u) |^. In NLO, the evolution 


of the wave function is given by 0 : 

/ 

V3^(M,/i) = ^a„(/ro)exp - / 

” \ Os(Mo) 



k>n 


(31) 


with oo = 1. The coefficients are due to mixing effects, induced by the fact that the 

polynomials 4 /„(m) are the eigenfunctions of the LO, but not of the NLO evolution kernel. 
The QGD beta-function (3 [O and the anomalous dimension y"’ of the n-th moment 


of the wave function have to be taken in NLO. Explicitly [|I5 

2 


a.. 




Att) 


71 


(32) 


with 


7o° 

7o 


7o 


0 , 

100 

728 

1 ^’ 


7i = 0 , 

34450 


7? = 


243 
_ 4 _ 662846 

“ 3375 


830 

- Nf, 

81 ’ 

31132 
~ 2025 


N, 


(33) 


The NLO mixing coefficients are P, 1^ 


<(a) = 


M 


nk 


7S - 7o" - 2A 


/ 

1 - 

/ asifi) 

7n"-7^-2/3o\ 
\ 2/30 

V 

^^^(/ro) 

) 

/ 


(34) 





































where the numerical values of the hrst few elements of the matrix Mnk are 


Mo2 = - 11.2 + 1.73Nf, Mo 4 = -1.41 + 0.565iVF, M 24 = - 22.0 + 1.65Nf. (35) 

With the above input and (|3TD we hnd a 2 (/iB) = 0.218 and = 0.084. 

Now, we are ready to exploit the sum rule (^) numerically. In Fig. 2, the prod¬ 
uct /b/"''( 0) is plotted as a function of the Borel parameter M^. The 0{as) correction 
turns out to be large, between 30% and 35% , and stable under variation of M^. More 
specihcally, in the interval = 8 -i -12 GeV^ we obtain in LO 


and in NLO 


/^/+( 0 ) = 0.0229^0.0224 

GeV, 

/^(O) 

= 0.163^0.160, 

(36) 

/^/+( 0 ) = 0.0306 ^ 0.0295 

GeV, 

r(0) 

= 0.170^0.164 

(37) 


where Jb = 140 MeV and 180 MeV has been used, respectively. Note the almost complete 
cancellation of the NLO correction in f~^. Furthermore, Fig. 3 shows the momentum 
dependence of the form factor in the region 0 < < 15 -h 17 GeV^ for = 10 

GeV^, where the sum rule (^) is expected to be valid. Finally, it is interesting to compare 
the fi dependence in LO and NLO. This is done in Fig. 4. The very mild fi -dependence 
in LO only results from the evolution of the wave function. In NLO, the p-dependence 
is stronger than in LO but similar to the /x-dependence of /b- As a result, the residual 
scale dependence of /■*■ is again mild. 

The above results refer to the leading twist 2 approximation. If one adds the LO twist 
3 and 4 contributions, one obtains at = 0 


/+(0) = 0.27 . 

This value should be compared with the LO estimate = 0.30 obtained in 


(38) 


5. In this paper, we have presented the perturbative QGD correction in 0{as) 
to the leading twist 2 approximation of the light-cone sum rule for the i? —> vr form 
factor f~^. Both UV and collinear divergences are handled by dimensional regularization 
and MS renormalization. The collinear divergences in the hard amplitude are factorized 
and absorbed in the evolution of the light-cone wave function. Numerically, the 0{as) 
correction to the product /s/"'" amounts to about 30%. We have shown that this large 
correction is almost completely compensated by the corresponding correction to the two- 
point sum rule for Jb- The remaining 0{as) effect on is therefore small. This hnding 
improves the accuracy and reliability of the light-cone sum rule estimate substantially. 

Furthermore, we have shown that the 0{as) correction to the sum rule for /b/"'" 
depends only very weakly on the momentum transfer. This observation, together with 
the above-mentioned compensation strongly suggests that the dominant 0{as) effect in 
the correlation function comes from the 75 vertex (see Fig. Ic) which is also present 
in the two-point correlation function for /b. Thus, our calculation strongly supports the 
conjecture i, I, pT| that the perturbative correction may drop out in the ratio /b/^//b. 

Recently, an estimate of the perturbative correction to the B ^ tt form factor was 
obtained m in a different approach combining the constituent quark model for B and tt 
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with light-cone wave functions. Although the results agree qualitatively it is difficult to 
directly compare our result with this model-dependent calculation. 

A more detailed account of our calculation as well as applications to various exclusive 
B and D decays will be published elsewhere. 

Acknowledgements. 

We are grateful to A. Ali, V. Braun, A. Grozin and A. Vainshtein for useful discussions. 
This work is supported by the German Federal Ministry for Research and Technology 
(BMBF) under contract number 05 7WZ91P (0). 


10 



References 


[1] V.L. Chernyak and A.R. Zhitnitsky, JETP Lett. 25 (1977) 510; Sov. J. Nucl. Phys. 
31 (1980) 544; Phys. Rep. 112 (1984) 173. 

[2] A.V. Efremov and A.V. Radyushkin, Phys. Lett. B94 (1980) 245; Teor. Mat. Fiz. 
42 (1980) 147. 

[3] G.P. Lepage and S.J. Brodsky, Phys. Lett. B87 (1979) 359; Phys. Rev. D22 (1980) 
2157. 

[4] V.M. Belyaev, A. Khodjamirian and R. Riickl, Z. Phys. C 60 (1993) 349 . 

[5] V.M. Belyaev, V.M. Braun, A. Khodjamirian and R. Riickl, Phys. Rev D 51 (1995) 
6177. 

[6] A. Khodjamirian and R. Riickl, in “Continuous Advances in QCD 1996”, 
edited by M.I. Polikarpov ( World Scientihc, Singapore, 1996), pp. 75-83; 

A. Khodjamirian, R. Riickl and Ch. Winhart , in preparation . 

[7] F. del Aguila and M.K. Chase, Nucl. Phys. B193 (1981) 517. 

[8] E. Braaten, Phys. Rev. D28 (1983) 524. 

[9] E.P. Kadantseva, S.V. Mikhailov and A.V. Radyushkin, Sov. J. Nucl. Phys. 44 
(1986) 326. 

[10] F.M. Dittes and A.V. Radyushkin, Phys. Lett. B134 (1984) 359; 

M.H. Sarmadi, Phys. Lett. B143 (1984) 471; 

S. V. Mikhailov and A.V. Radyushkin, Nucl. Phys. B254 (1985) 89. 

[11] D.J. Broudhurst and S.C. Generalis, preprint OUT-4102-8/R (1981) (unpublished); 
D.J. Broadhurst, Phys. Lett. BlOl (1981) 423; 

T. M. Aliev and V.L Eletsky, Sov. J. Nucl. Phys. 38 (1983) 936. 

[12] Particle Data Group, Phys. Rev. D54 (1996) 1. 

[13] V.M. Braun and I.E. Filyanov, Z. Phys. C48 (1990) 239. 

[14] V.M. Braun and I.E. Filyanov, Z. Phys. C44 (1989) 157. 

[15] A. Gonzales-Arroyo, C. Lopez and F.J. Yndurain, Nucl. Phys. B153 (1979) 161. 

[16] P. Ball, V.M. Braun and H.G. Dosch, Phys. Lett. B273 (1991) 316. 

[17] A. Szczepaniak, Phys. Rev. D54 (1996) 1167. 


11 










Figure 1: Feynman diagrams contributing to the correlation function (j^): (a) zeroth order 
in tts, (b-g) hrst order in as- 
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Figure 2: Light-cone sum rule estimate for /b/^( 0) in leading twist 2 approximation as 
a function of the Borel parameter NLO (solid ) in comparison to LO (dashed). 



Figure 3: Momentum dependence of the form factor /^(p^) in leading twist 2 approxima¬ 
tion; LO (dashed) in comparison to NLO (solid). 
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Figure 4: Scale dependence of the light-cone sum rule estimate of /b/’''( 0) in leading twist 
2 approximation: NLO (solid) in comparison to LO (dotted). 
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